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For hole systems with an effective spin j — 3/2, we present an invariant decomposition of the spin 
density matrix that can be interpreted as a multipole expansion. The charge density corresponds 
to the monopole moment and the spin polarization due to a magnetic field corresponds to a dipole 
moment while heavy hole-light hole splitting can be interpreted as a quadrupole moment. For quasi 
two-dimensional hole systems in the presence of an in-plane magnetic field B\\ the spin polarization 
is a higher-order effect that is typically much smaller than one even if the minority spin subband is 
completely depopulated. On the other hand, the field Bu can induce a substantial octupole moment 
which is a unique feature of j = 3/2 hole systems. 

PACS numbers: 73.21. Fg, 72.25.Dc, 71.15.Mb 



c3 



i 

-a 
c 

o 
o 



> 
o 

m 

(N 
O 

o 

-I— > 



i 

C 

o 
o 



X 



I. INTRODUCTION 

A powerful approach to the quantum many-body prob- 
lem is density functional theory. Hohenberg and Kohn 1 
showed that the (nondegenerate) ground state of a sys- 
tem of interacting particles is a unique function of the 
corresponding ground state density p. The ground state 
of the system is thus fully characterized by its density p 
(Ref. 0). For electrons with spin j = 1/2 the density p 
must be replaced by a 2 x 2 spin density matrix p. It is 
well-known that p can be parameterized by the density 
p = tr p and (the magnitude of) the spin polarization. 2 
The spin polarization at a magnetic field B > is thus 
an important quantity for our understanding of many- 
particle phenomena. Recently, the interest in this sub- 
ject has been renewed because experiments have indi- 
cated that the spin polarization of quasi two-dimensional 
(2D) systems due to a magnetic field Bu parallel to the 
2D plane affects the apparent metallic behavior of these 
systems. 3 Furthermore, the spin polarization is an im- 
portant parameter for possible applications in the field of 
spintronics* In quasi two-dimensional (2D) systems the 
spin polarization has been studied by applying a mag- 
netic field B\\ parallel to the 2D planei 5 i 6 i 7 i 8 i 9 i 10 i n i 12 i 13 i 14 
If the Zeeman splitting due to -By becomes sufficiently 
large, the minority spin subband is completely depop- 
ulated. In such a situation, electron systems are fully 
spin-polarized. This approach has been used by several 
groups to study the spin polarization in low-density 2D 
electron systems LS*2iiS 

Recently, there has been considerable interest in the 
spin polarization of quasi 2D hole systemsii*i2iiii4. Hole 
systems in the uppermost valence band of many common 
semiconductors such as Ge and GaAs are different from 
electron systems due to the fact that the hole states have 
an effective spin j = 3/2 (Ref. ITlj) . Usually, hole states 
in quasi 2D systems are discussed using a basis of angu- 
lar momentum j = 3/2 eigenfunctions with quantization 
axis perpendicular to the 2D plane. Subband quantiza- 
tion yields so-called heavy hole (HH) states with z com- 
ponent of angular momentum m = ±3/2 and light hole 
(LH) states with m = ±1/2. Often it is assumed that 



the spin degree of freedom of these HH and LH states 
behaves similar to the spin of j — 1/2 electron states. 
However, recently it has been shown that the spin po- 
larization of quasi 2D HH systems in the presence of an 
in-plane magnetic field behaves rather different from the 
more familiar case of j = 1/2 electron systems^ For ex- 
ample, the spin polarization of quasi 2D HH systems can 
change its sign at a finite value of B» . 

Similar to the case of j = 1/2 electron systems, the 
fundamental object for the characterization of j = 3/2 
hole systems is the spin density matrix p from which 
we can derive, e.g., the spin polarization. The domi- 
nant character of the occupied eigenstates of hole sys- 
tems depends on the quantization axis of the underlying 
j = 3/2 basis functions. Obviously, all observable quanti- 
ties such as the spin polarization may not depend on this 
choice. Therefore, it is necessary to formulate the spin 
density matrix in a way such that observable quantities 
can be calculated independent of the particular choice 
for the basis functions that are used. It is the goal of 
this paper to present such a theory. Using the theory of 
invariantsi2ii& we will derive an invariant decomposition 
of the spin density matrix of j = 3/2 hole systems that 
can be interpreted as a multipole expansion. 



II. INVARIANT DECOMPOSITION OF THE 
SPIN DENSITY MATRIX 

A. General theory 

In general, the single-particle states in quasi 2D sys- 
tems with effective spin j (the Kohn-Sham states) can be 
characterized by multicomponent envelope functions^ 
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(1) 



Here £™ a t.n ( z ) denotes the mth component of the enve- 
lope function that modulates the band-edge Bloch func- 
tion u m {r), m labels the z components of j, a is the 
subband index, a is the spin index, and fen = (k x ,k y ,0) 
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is the in-plane wave vector. The (2j + l)-dimensional 
Hermitian, positive definite spin density matrix p has 
then the matrix elements (temperature T = 0) 

(m|m')= [ dk f\ 0[±(£F-£aa(fc||))] 

Q.fJ 

x /^C*„ (*)££*,,(*)■ (2) 

Here E'ao^fcn) is the spin-split subband dispersion and 
i?F is the Fermi energy. The upper (lower) sign in Eq. 
© applies to electron (hole) systems. We integrate over 
all orbital degrees of freedom that are irrelevant in the 
present discussion. The density is given by the trace of 
p. Obviously, the trace is independent of the particular 
choice for the set of basis function u m (r). 

For 3D homogeneous hole gases {B — 0) the matrix 
p is proportional to the unit matrix, i.e., we have equal 
contributions of HH and LH states because we have no 
preferential direction for the quantization axis of angu- 
lar momentum^ (This was overlooked in Ref. l2fT.) For 
inhomogenous systems such as quasi-2D systems, the ma- 
trix p has a more complicated structure, even at B = 0. 
In principle p contains all information about the ground 
state of the system. However, as the Hermitian 4x4 ma- 
trix p is characterized by 16 independent parameters, it 
appears impractical to characterize a j = 3/2 system by 
its density matrix. Thus, Bobbert et al£L suggested to 
parameterize the mixed HH-LH character of inhomoge- 
nous hole gases by the ratio between the LH and HH 
Fermi wave vectors. Enderlein et alw^ proposed to use 
partial densities for HH and LH states. The latter ap- 
proach was criticized by Bobbert et alm& because the 
partial hole densities depend on the set of basis func- 
tion u m (r). Recently, Karkkalnen et al~ suggested to 
treat mixed HH-LH systems as two-component systems 
where the HH and LH components are essentially inde- 
pendent of each other. Until today, the interpretation of 
the 4x4 spin density matrix p of particles with j = 3/2 
has remained an open problem. 

We show here that the theory of invariants' 7 *^ allows 
to decompose p for systems with arbitrary point group 
Q into symmetry-adapted irreducible tensors which pro- 
vide a rigorous yet physically transparent classification of 
the different terms in p — independent of the particular 
choice for the set of basis functions. We assume that the 
system is characterized by an N x N multiband Hamilto- 
nian Ti. The basis functions of TL transform according to 
an iV-dimensional representation T a of Q. For simplic- 
ity, we assume that the representation T a is irreducible. 
The density matrix p of this system can be written in 
the form 

P = . (3) 

Here are N x N matrices with all elements equal 
to zero, except for the element {p 1 v) that equals one. 
The components {p^} and the matrices {M M „} trans- 
form according to the product representation T a x T* . 



In general, this representation is reducible, i.e., we have 

r a xr*=^i>. (4) 



Equation implies that we can decompose the set of 
matrices {M^} (the components {p^v}) into irreducible 
tensors {Mp} ({pp}) that transform according to the 
representations Tp of Q . The symmetry-adapted expan- 
sion of the density matrix p is then given by 17 

p = Y J M t3 -pp=Y J M ^Ph^ ( 5 ) 

/3,A 

where A labels the elements of the irreducible represen- 
tation Tp. The terms pp\ can be interpreted as the pro- 
jection of p on the basis matrices Mp\. If the matrices 
{Mp\} are orthonormalized we thus have 

P/3A = tr(M /3A p) ■ (6) 

According to Eq. JSJ), we get a decomposition of the 
density matrix p in terms of invariants Mp ■ pp. Of- 
ten, these invariants provide already a clear interpreta- 
tion of the different terms in p. A second advantage 
of the invariant expansion (0 is the following. Using 
the Clebsch-Gordan coefficients for the irreducible repre- 
sentations {r^} of Q we can construct scalar quantities 
{p s } from the irreducible tensors {pp} which transform 
according to the unit representation Fi (Ref. Il7l) . These 
scalars are sums of products of the elements {pp\}- By 
definition, they are independent of the particular choice 
for the set of basis functions that was used to evaluate 
the spin density matrix. Therefore, the scalars {p s } rep- 
resent a convenient set of independent variables for a pa- 
rameterization and characterization of the spin density 
matrix of a system with point group Q . 

B. Multipole expansion of the density matrix for 
systems with point group SU(2) 

For a system with spherical symmetry we have Q = 
SU(2). We denote the irreducible representations of 
SU(2) by 2?j where j is the angular momentum. Here 
the invariant expansion (jSJ) is equivalent to a multipole 
expansion of the density matrix. It follows from the well- 
known relation for the addition of angular momenta^ 

)+)' 
)"=|j-j'l 

that for a system with angular momentum j the invariant 
expansion (JSJ) contains multipole moments p } up to the 
order 2j. The squared magnitude of the jth multipole 
moment is given by2» 

pf=Yl (-l) m /Wj,-m. (8) 
m=-j 
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TABLE I: Invariant decomposition of the spin density matrix 
of a system with spin j = 1/2. 
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75 12x2 



[a x + ia y ) 



75 (( 2 I 2 / \ 2 I 2; 
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TABLE II: Tensor operators /Cj for the point group SU(2). 



Vj m /Cj r 



£>o fc 2 = {fc + ,fc_} + fc 2 

B 2 
V 2 2 -ifc?. 

1 -{fc+,fc z } 

-J=(2fc 2 -{fc +J fc_}) 



T> 3 3 



2 ^(fc 2 ^ + 2{fc+,fc z }B + ) 

1 ^=(4 B- + 2{k+, k-}B+ - 4k 2 z B+ - 8{k+, k z }B z 
^[(2k 2 z - {k + ,k.})B z - {k+,k z }B. - {k.,k z }B, 



By definition, these scalars are independent of the quan- 
tization axis of angular momentum that is used for the 
basis functions u m {r) in Eq. 



1. j ~ 1/2 Electron systems 

First we apply the above formalism to the well-known 
case of a density matrix for systems with angular mo- 
mentum j = 1/2. Using standard angular momentum 
algebra^ we obtain the invariant decomposition listed 
in Table [IJ where a x , a y , and a z are the Pauli spin 
matrices. We have omitted the terms with m < 
which can be obtained using Mj._ m = (— l)' ,+m Mj m and 
Pj- m = (-l)> +m p* m (Ref.il). We note that the terms 
Pj m can be obtained either from Eq. 10 or by using the 
Clebsch-Gordan coefficients of SU(2) (Ref-H). By sum- 



ming over the invariants Mj ■ pj it is easy to check that 
these invariants provide a decomposition of the spin den- 
sity matrix p. 

Using Eq. (jHJ we obtain the scalars 



Pi 



v/2 



(9a) 



M<^M4l^» 2 + 2|U|-±)| 2 .(9b) 



Apart from a prefactor 1/V2, the magnitude of the 
monopole moment, is equal to the density. (We re- 
mark that by definition po is already a positive scalar 
that transforms according to T>q so that we always have 
Po = \fPo = Po-) Apart from a prefactor 1/2, the 
(squared) magnitude of the dipole moment, p\ , is equal 
to the (squared) magnitude of the polarization. Due to 
the rotational invariance implied by Q = SU(2), the den- 
sity matrix is completely characterized by po and p\ . 

Consistent with the invariant expansion of the density 
matrix one can perform an invariant expansion of the 
Hamiltonian 7^2x2 of j = 1/2 electron systems. In lowest 
order of k and B = (ih/e)k x k we obtain the tensor 
operators listed in Table ILTI Here k± = k x ± ik y , and 
{. . .} denotes the symmetrized product of its arguments, 
e.g., {A, B} = \{AB + BA). We then obtain 



"^2x2 — 



^2r. 



-Mo-7Co+^=MBM r /Ci+U(r)l 2x2 . (10) 
v2 



Here V(r) is the external potential. The kinetic en- 
ergy operator with effective mass m* transforms like a 
monopole and the Zeeman term with effective g factor 
g* transforms like a dipole. The symbol pb denotes the 
Bohr magneton. We see here that for j = 1/2 systems 
the jth invariant in characterizes the dynamics of 

the jth invariant of the spin density matrix. However, 
we want to emphasize that this simple correspondence 
between the different terms in the Hamiltonian (|10fl and 
the multipole moments © is valid only because we can 
factorize the eigenfunctions of 7i 2 x2 into an orbital part 
and a spin part. In general, no such factorization can 
be performed for j — 3/2 hole states. This is the reason 
why we will find below that hole systems can have an 
octupole moment even though the corresponding Hamil- 
tonian does not contain any octupole term. 



2. j — 3/2 Hole systems 



We now perform the invariant decomposition of the 
spin density matrix for a system with j — 3/2. According 
to Eq. J7|) it contains multipole terms up to j = 3. Using 
standard angular momentum algebra™ we obtain the in- 



variant decomposition listed in Table ITTT1 Here J x , J, 



and J z are the angular momentum matrices for j = 3/2 
and J± = J x ± iJ y . Once again, the (squared) magni- 
tudes of the multipole moments can readily be calculated 
by means of Eq. JHJ . We obtain 
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+ I(f|-|)-(I|-I)| 2 + 2| 


2I-2M- 
(ll-l)l 2 - 


3 1 3^ 
2 1 2 / 
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Apart from a prefactor 1/2, the monopole moment pq is 
equal to the density, and apart from a prefactor 9/20, the 
squared dipole moment p 2 is equal to the squared mag- 
nitude of the polarization. 25 The quadrupole moment 
Pi and the octupole moment /9 3 have no equivalent in 
j = 1/2 systems. The magnitude P2 of the quadrupole 
moment quantifies the effect of HH-LH splitting. If only 
the diagonal elements (m\m) of p are nonzero, pi is equal 
to the difference between the partial HH and LH densities 
used by Enderlein et alm^ (apart from a prefactor 1/2). 
For j — 3/2 systems, a complete characterization of the 
density matrix by means of scalar quantities requires to 
take into account additional terms beyond those obtained 
from Eq. (JSJ. They characterize, e.g., the interaction be- 
tween different multipole moments. It follows from Eq. 
(|7|) that these scalars are of higher order in (m\m'). 

The dynamics of j = 3/2 hole systems is characterized 
by the 4x4 Luttinger Hamiltonian 7^4x4 fRef. Il5^ . Using 
Tables ITT1 and ITTT1 we express 7^4x4 in terms of spherical 
tensor operators 2 ^! 



TLa 



h 2 

—M • TCo - 2VbKp B Mx ■ TCi 
m 



| ~ H 2 2 72 + 3 73 

mo 5 
+U(r)l 4 x4 , 



Mi ■ K 2 + CM 3 • K 3 



(12) 



where 71, 72, and 73 are the Luttinger parameters and 
k is the isotropic effective g factor. We neglect here the 
small terms with cubic symmetry. They are considered in 
the numerical calculations discussed below. Once again, 
the magnetic field is acting in spin space as a dipole field 
(the operator TCi). Obviously, the quadrupole and the 
octupole term in 7^4x4 have no equivalent in j = 1/2 
systems. The quadrupole operator K-i is responsible for 
the HH-LH splitting^ We note that in bulk material, 
uniaxial or biaxial strain also gives rise to a splitting of 
HH and LH states because analogous to 1Ci the strain 
tensor contains a term 7C 2 that transforms accordi ng t o 
T>2 so that we get another invariant cx Mi ■ K! 2 (Ref. |27j). 
In Eq. (|12fl . the symbol £ denotes the prefactor of the 
octupole term. Starting from the 8x8 Kane model we 
obtain in fourth order perturbation theory^ 



V2 e P 4 
3 



' - h El 1 - 



1 

£0 



1 

A 



(13) 



where P is Kane's momentum matrix element, E is the 
fundamental gap, and Ao is the spin-orbit split-off gap. 
If we insert typical values for P, Eq, and Ao it can be 
seen that £ is rather small so that the octupole term in 
7^4x4 can safely be neglected. However, we show below 
that, nonetheless, the spin density matrix of hole systems 
can have a large octupole moment. 



III. SPIN DENSITY MATRIX OF QUASI 2D 
SYSTEMS 



A. j — 1/2 Electron systems 

As an example we consider next the spin density ma- 
trix of quasi 2D electron and HH systems in the presence 
of an in-plane magnetic field Bi\. For j = 1/2 electron 
systems at Bn > we can choose the eigenstates inde- 
pendent of Bu 



2- 



-( 1 

V2 U 1 



(14) 



Using Eq. (|9b|) we then obtain for the dipole moment 
normalized with respect to the total 2D density 



1 / m*g*p B B\\ \- 
- ; 2 I. ^Ir\ ) 



Bu < B D 



B n > Br, 



(15) 



where we have assumed that only the lowest orbital sub- 
band a = 1 is occupied, N denotes the total density 
in the 2D system, and B p = 2irh 2 N / (m* g* p B ) is the 
magnetic field at which the system becomes fully spin- 
polarizedi2& In Eq. i|14|) the index z at the spinor indicates 
that we have chosen the spin quantization axis in z direc- 
tion perpendicular to Bu . Of course, the dipole moment 
Pi may not depend on this seemingly inconvenient choice. 
It is easy to check that we get the same result if the spin 
quantization axis is chosen parallel to B11 . 
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TABLE III: Invariant decomposition of the spin density matrix of a system with j — 3/2. 
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B. j — 3/2 Heavy hole systems 

Next we consider a 2D HH system in the presence of 
an in-plane magnetic field Bn . We write the degenerate 
HH eigenstates at B = in the form 16 



K> = 



2tt 



V2 



/ 1 \ 




V±i/ 



(16) 



Here we have assumed that the Luttinger Hamiltonian 
(|12[l is expressed in a basis of j = 3/2 angular momen- 
tum eigenfunctions in the order m = 3/2, 1/2, —1/2, and 
—3/2; and the index z at the spinor indicates that the 
quantization axis of these basis functions is perpendic- 
ular to the 2D plane. For simplicity we have neglected 
HH-LH mixing due to a nonzero in-plane wave vector fen 
and we also neglected the fen dependence of the envelope 
functions £^ H (z). These aspects are fully taken into ac- 
count in the numerical calculations discussed below. A 
finite in-plane magnetic field Bn gives rise to a mixing of 
HH and LH states. In first order degenerate perturba- 
tion theory we obtain the following expressions^ for the 
quasi 2D HH states in an in-plane field Bn 



where 



K 
Q 



2tt Ka [ ' V2 



K/lB-Bll 



±1 



. (17) 



E h — E l ' 



2 72 + 3 73 2mo<£?Vl£« H } (A* B B||) 



h 2 



E h - E l 



(18a) 



(18b) 



Herei£ (E l a ) is the energy of the ath HH (LH) subband 
at B11 = 0, and we have used the gauge A = (0, — zBn, 0). 
The expression for Q in Ref . differs slightly from our 
Eq. i|18b(l because we have used here the spherical ap- 
proximation (fT2*|l of H4x4- As discussed in Ref.Hfilwe can 
often treat the HH-LH mixing characterized by JC and Q 
as a small parameter even if Bn is so large that only the 
majority spin subband is occupied. This is a consequence 
of the frozen angular momentum of HH states in quasi 
2D systems. 

Using Eq. I jlljl we then obtain for the multipolc mo- 
ments normalized with respect to the total 2D density 



pi 



Pi 



5 



1 



3/C 2 [(MO 2 - 1] 



(SN) 2 



K, 1 



3(<W) 5 



1 8 

— QK SN 

5 



9 



g 2 (sn) 1 



(19a) 
(19b) 

r 

(19c) 



a ot 



where SN = (N + — N-)/(N + + N-) denotes the normal- 
ized difference between the 2D densities N± in the major- 
ity and the minority spin subband (with N = N + + N— ) , 
and we have assumed that only the lowest HH subband 
q = 1 is occupied (which is usually the case for 2D hole 
systems). We note that prior to deriving the multipolc 
moments <|19|1 one must normalize the perturbed wave 
functions (|17fl . 

HH-LH mixing gives rise to a dipole moment p\ pro- 
portional to the small parameters K, and Q (Ref. Hq) . 
In lowest (i.e., zeroth) order of the HH-LH mixing, the 
quadrupole moment p 2 , is equal to 1/4, i.e., p 2 is es- 
sentially independent of Bn, which is a consequence of 
HH-LH splitting. Finally, we also get a substantial oc- 
tupole moment even though we assumed Q = in the 
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Luttinger Hamiltonian 1|12|) . If the minority spin sub- 
band is completely depopulated the normalized octupole 
moment equals 1/2 in lowest (i.e., zeroth) order of the 
HH-LH mixing. We note that the Zeeman energy split- 
ting of the HH states in the presence of an in-plane mag- 
netic field -By is proportional to B^ fR.ef.ll8l). Neglecting 
the nonparabolicity of the HH energy dispersion we thus 
have 5N oc Bj. 

It is interesting to trace back the physical origin of 
the multipole moments P2 and p^. In spite of the fact 
that we do not have a simple correspondence between the 
multipolcs in the Luttinger Hamiltonian [Eq. (|12[) ] and 
the multipole moments pi of the spin density matrix [Eq. 
(|19J) ] , we get nonzero multipole moments pi and P3 only 
if the Hamiltonian (|12fl contains a quadrupole term oc 
M2 ■ JC2 ■ If the quadrupole term in Eq. I|12f> were zero the 
HH eigenstates (f lfif) at B — were degenerate with the 
corresponding LH states. Similar to the electron states 
(fH|l these four degenerate states could be factorized into 
a scalar orbital part times a spinor that depends only on 
B, and we would have P2 = pz = 0. To the best of our 
knowledge, the quadrupole term cx M2 ■ K-2 is a unique 
feature of the Hamiltonian for the extended valence band 
states in semiconductors 4^SSiSl 

It is instructive to check that the scalars pi in Eq. 
(|19f) are indeed independent of the quantization axis e of 
angular momentum used in Eqs. l|16fl and l|17fl . If e is 
chosen parallel to Bi\ Eq. I|17|l must be replaced by22i 



ik\\-r\\ 



I a+ / 



l*™> 



/ -1 + 3JC + 3G \ 


V3(l + /C + g) 




(20a) 



1 r ° ^ 



\ 1 + 3/C - 3G / 



20b) 



In this basis the dominant spinor components of the HH 
states have \m\ — 1/2. Nonetheless, we obtain the mul- 
tipole moments (I19|) . 



C. Numerical results 

To confirm the qualitative results in Eq. I|19|l we show 
in Fig. ^ the self-consistently calculated multipole mo- 
ments p\ as a function of the in-plane magnetic field -By 
for a symmetric (100) GaAs-Alo.3Gao.7As quantum well 
with hole density N = 5 x 10 10 cm~ 2 and well width 
w = 150 A. The numerical calculations follow Ref. 0. 
We evaluate Eq. @ by means of analytic quadratic Bril- 
louin zone integration^ The solid lines in Fig. 2] have 
been obtained using the spherical approximation (|12fl of 
the Luttinger Hamiltonian 7^4x4- However, we have ne- 
glected the small octupole term in Eq. I|12f> . For com- 
parison we also show the results based on the 8x8 Kane 
Hamiltonian 7^8x8 which includes implicitly the octupole 
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FIG. 1: (a) Normalized spin subband density N-/N of 
the HH minority spin subband and (b-d) normalized mul- 
tipole moments p\ as a function of the in-plane magnetic 
field B\\ calculated self-consistently for a symmetric (100) 
GaAs-Alo.3Gao.7As quantum well with hole density iV = 
5 x 10 10 cm -2 and well width w = 150 A. The solid lines 
have been obtained by means of the spherical approximation 
(1 1 21) of the Luttinger Hamiltonian. The dashed (dotted) lines 
correspond to the spherical (cubic) approximation of the 8x8 
Kane Hamiltonian, respectively. 



term. The dashed lines in Fig. ^ are based on the spher- 
ical approximation of Tisxs while the dotted lines also 
take into account the cubic terms^ 

We can see in Fig. ^ that the complete depopulation 
of the HH minority spin subband due to By does not im- 
ply full spin polarization of the system*^ In Fig. (solid 
lines), the minority spin subband is completely depopu- 
lated at B D w 23.9 T while p\{B D ) « 0.0093. The latter 
value is much smaller than 9/20 = 0.45, the value of p\ in 
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a fully spin-polarized HH system. The zero of p\(B\\) at 
B« ~ 18.4 T reflects a sign reversal of the spin polariza- 
tion vector which is a unique feature of 2D HH systems. 16 
The derivative of pf(Bn) is discontinuous at Bu = Bp- 
At -By = the quadrupole moment p\ is slightly smaller 
than 1/4, the value of p\ in a pure HH system. This is 
a consequence of the fey -induced HH-LH mixing which 
was fully taken into account in Fig. ^ For Bu > we 
observe only a small decrease of p\. This is due to the 
fact that the HH-LH splitting E\ -E\k, 6.7 meV [i.e., 
the denominators in Eq. (|18|) ] is the largest energy scale 
in the system so that the HH states have a frozen angular 
momentum perpendicular to the 2D plane. For compari- 
son, we note that the Zeeman energy splitting of the HH 
subband at Bu is ~ 0.4 meV. The spin polarization due 
to Bu is thus a higher-order effect in 2D HH systems. 

It is remarkable that the octupole moment p§ at By = 
Bo is close to 1/2, the largest possible value of p| in 
a 2D HH system. This value is essentially independent 
of whether we use 7^4x4 without an octupole term or 
whether we use 7isx8 which includes implicitly the oc- 
tupole term oc M% ■ K3. For Bu > Bp, the octupole 
moment p^{Bu) remains essentially constant, consistent 
with Eq. (|19c|l . These findings suggest that an in-plane 
magnetic field -By provides an efficient tool to study 2D 
HH systems with a large octupole moment but with a 
small dipole moment (i.e., with a small spin polariza- 
tion). We can obtain a 2D HH system with a large spin 
polarization but with a small octupole moment if the 
magnetic field is applied perpendicular to the 2D plane. 

Figure ^ shows that the cubic terms are only small 
corrections in Hsxs- They can safely be neglected in a 
discussion of the main features in Fig. ^ This is con- 
sistent with the fact that the cubic terms in 7isx8 (and 
Ti.4x4) are proportional to 73 — 72 which is typically a 
small parameter^ On the other hand, it was important 
for the results shown in Fig.^tbat the numerical calcula- 
tions took into account HH-LH mixing due to a nonzero 
in-plane wave vector fey . As the analytical expressions in 
Eq. Q19(l neglect this effect, we do not present a direct 
comparison between Eq. I|19fl and the numerical results 
in Fig. ^ We remark that it is straightforward to include 
the fey -induced HH-LH mixing in Eq. i|19|) . However, we 
do not reproduce here the lengthy formulas. 



IV. OTHER SYSTEMS AND OUTLOOK 

The numerical calculations shown in Fig. ^ have ne- 
glected many-particle effects beyond the Hartree approx- 
imation. These effects will be important for a quantita- 
tive analysis of the spin multipole moments in the density 
regime discussed here. However, to the best of our knowl- 
edge an appropriate theory for particles with j = 3/2 is 
currently not available. It represents a formidable task to 
develop such a theory. We note that the symmetry argu- 
ments developed here are not affected by many-particle 



effects. The multipole expansion of the spin density ma- 
trix can thus provide a useful starting point for a more 
quantitative theory. We expect that the essential features 
in Fig. [2 will be confirmed by such an investigation. 

Our results can readily be generalized to spherically 
symmetric systems with arbitrary j, both half integer 
and integer j24 as well as to systems characterized by the 
crystallographic point groups— As an example, we con- 
sider a j = 1/2 electron system in a spatial environment 
with point group Q = C^v Here the 2x2 spin density 
matrix is completely characterized by five independent 



scalars 

po = + ( 21a ) 

P 1 = IGI-!> + (-!li>l 2 (21b) 

p 2 = \(hH)-(-m 2 (2ic) 

^ = i<iii>-<-ii-i>r (2W) 

p 4 = «ili>-(-||-l»«l|-l> + <-||l» 

*«*!-*> -<-*!*»■ (2ie) 



Each term in Eq. (|21|) has a clear physical interpreta- 
tion. Once again, the scalar po is the charge density. The 
scalars p\ , pi , and p 3 describe the magnitude of the spin 
polarization in x, y, and z direction, respectively. These 
quantities must be distinguished from each other in a 
system with point group Ci v . Finally, P4 describes the 
relative sign of the spin polarization in x, y, and z direc- 
tion. The five scalars pi thus form a decomposition of the 
four independent parameters of a Hermitian 2x2 matrix 
into their magnitudes and their relative signs. For the 
explicit formulas in Eq. (12 1H we have assumed that the 
angular momentum quantization axis of the basis func- 
tions is parallel to the symmetry axis of the point group 
C2V ■ We emphasize that the magnitude of the scalars pi 
in Eq. (|2"TJ) does not depend on this choice. 

Finally, we want to briefly address the question of how 
one can measure the multipole moments of the spin den- 
sity matrix of hole systems. The polarization of j = 1/2 
electron systems can be probed by measuring the hyper- 
fine interaction between the electrons and the atomic nu- 
clei. Electrons in the conduction band of common semi- 
conductors like GaAs originate in s-like atomic orbitals. 
Therefore, the electrons have a large probability density 
at the atomic nuclei so that the hypcrfinc interaction is 
large. The j = 3/2 hole states in the valence band, on 
the other hand, originate in p-like atomic orbitals which 
have a vanishing probability density at the atomic nuclei. 
The hyperfine interaction is therefore much weaker— so 
that it will be difficult to use this technique to measure 
the polarization of j = 3/2 hole states. 

Indirect evidence for the unusual features of the spin 
density matrix of hole systems can be obtained by varying 
the confinement of quasi 2D HH systems perpendicular 
to the 2D plane. This will change the HH-LH splitting 
E% — E l a in the denominators of Eq. l|18f) . Here our re- 
sults are in good qualitative agreement with the trends 
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observed in recent experimentsiii*i2ii& A detailed com- 
parison with experimental results will be the subject of 
a future publication. 
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